Brane curvature and supernovae la observations 
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It is well known that modifications to the Friedmann equation on a warped brane in an anti de 
Sitter bulk do not provide any low energy distinguishing feature from standard cosmology. However, 
addition of a brane curvature scalar in the action produces effects which can serve as a distinctive 
feature of brane world scenarios and can be tested with observations. The fitting of such a model 
with supernovae la data (including SN 1997ff at z ~ 1.7) comes out very well and predicts an 
accelerating universe. 
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In recent times there has been considerable activity in 
large extra dimensions and the recovery of the four di- 
mensional General Relativity (GR) as an effective theory 
from a more fundamental theory. The basic idea in these 
theories is the existence of a higher dimensional bulk in 
which our universe, called 3-brane is sitting as a hyper- 
surface. The standard model matter fields are confined 
to the brane while gravity can, by its universal character, 
propagate in all (including extra) dimensions. The large 
extra dimensions also promise to solve the mass hierar- 
chy problem of the standard model of particle physics. 
A particular example is the model of Randall and Sun- 
drum jj, which has a warped extra dimension and has 
attracted a lot of interest. 

In this model the authors considered a 5-D anti de Sit- 
ter (AdS) bulk with a ^-symmetric extra dimension and 
by a proper fine tuning of bulk parameters they obtained 
a flat static brane with a vanishing effective cosmological 
constant on the brane. The standard Newtonian poten- 
tial is recovered in the model with r~ 3 corrections aris- 
ing from the massive Kaluza-Klein (KK) modes on the 
brane. The ^-symmetry is motivated by the reduction 
of M theory to E$ x Eg, heterotic string theory . 

There have been many generalizations of this model 
and various cosmological issues have been addressed 
therein ( see for example @). The recovery of Newto- 
nian force law in this kind of models is a non-trivial task 
and requires solving the perturbations of the bulk metric. 
It is not clear a priori whether Newtonian gravity can al- 
ways be recovered in such models. A counter example is 
provided by the conformally non-flat Nariai metric, for 
which there exists no massless graviton on the brane Q . 
On the other hand, this recovery has been established 
for the Schwarzschild-AdS (S-AdS) bulk (which is also 
conformally non-flat) with an FRW brane metric and it 
gives rise to various possible physical universes including 
the one with effective positive cosmological constant |5j. 
Hence it could be envisioned that our FRW universe may 
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be expanding out of a Schwarzschild black hole in an AdS 
bulk space. 

Using the Israel junction conditions [(J, one can write 
the analogue of the Friedmann equation on the brane 
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where G and A are respectively 4-D gravitational and 
cosmological constants, p is the matter density on the 
brane, S is the scale factor, a is the brane tension and 
C is the mass parameter of the bulk black hole. The 
four dimensional effective constants are related to the five 
dimensional bulk constants through the Israel junction 
conditions. 
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where A5 = — 6/Z 2 and G5 are the 5-D cosmological and 
gravitational constants. Here I is the radius of curvature 
of the bulk spacetime. 

The model differs from the standard cosmology in the 
following two terms in eq.Q). First, the density-squared 
term representing the local effects, which arises due to 
corrections in stress tensor by imposing junction condi- 
tions. This term decays as 1/S 8 in the radiation domi- 
nated epoch and would be insignificant even at the time 
of big bang nucleosynthesis (BBN) [3 • Second, the term 
varying as 1/S , commonly known as the dark radiation, 
represents the non local bulk effect in terms of mass of 
the black hole. This term enters into eq.Q) through the 
projection of bulk Weyl curvature on the brane and be- 
haves like radiation. It can affect the BBN and hence can 
be constrained by observations. As we shall show in the 
following, the dark radiation term can be safely neglected 
at the epoch even as remote in the past as z ~ 1.7, which 
is the highest redshift observed so far in the case of the 
supernovae la. Thus there survives no observational ef- 
fect of these brane world modifications and the model 
reduces to the standard FRW model without having any 
distinguishing features of its own. 

It is interesting to note that if one considers the brane 
models with Minkowski bulk 0, 0, the deviations 
from standard cosmology are expected at a cross over 
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scale r c defined as the ratio M|/2Mg, where M4 and M5 
are respectively the 4-D and 5-D Planck masses. The ac- 
tion in these models contains a term proportional to the 
curvature scalar on the brane which also helps in recovery 
of 4-D Newtonian gravity on the brane. The cosmologi- 
cal consequences of such an introduction in these models 
were first studied in Ref. [ToL ITU ] and the Friedmann 
equation with the cosmological constant on the brane in 
these models has also been derived [lJJ, EH ■ The compar- 
ison of these models with cosmological observations has 
also been pursued [l3j . 

The term proportional to curvature scalar can be 
viewed as a quantum correction and is also usually 
needed for a proper definition of stress-energy tensor on 
the boundary of S-AdS spacetime A variant of the 
AdS bulk brane model was considered by Kim, Lee and 
Myung |l5j in which they incorporated a brane curvature 
scalar term in the bulk-brane action via a small coupling 
parameter. This model provides a departure from the 
standard cosmology at late time evolution of the uni- 
verse. Here, our main aim is to test the model against 
the supernovae la observations and study its cosmologi- 
cal consequences. In fact, the model fits the SN la data 
very well. 

The brane in this model is a hypersurface given by the 
metric 
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embedded in a bulk composed of two patches of negative 
A5 spacetime whose action is 
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where K is the trace of the extrinsic curvature of the 
brane. The first term in the action is the standard 
Einstein-Hilbcrt action, the second term is the Gibbons- 
Hawking term which is necessary for the variational prob- 
lem to be well defined and the third term is the contri- 
bution of brane tension to the action. We are interested 
in looking at the incorporation of the brane curvature 
scalar term in the action coupling through the parameter 
b which is small so that its higher order terms in effective 
Einstein's equations on the brane can be neglected at the 
present e poc h. For this model the modified Friedmann 
equation [15J , in low energy limit (p <C er) , becomes 
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There would also be higher order terms in S in eq.JSJ, 
which have not been included for their effect being in- 
significant at the epochs of redshifts upto 1.7. The term 
a/S 2 would be dominant over the density and dark en- 
ergy terms at late times of expansion of the universe. It 
may be noted that 1/S 2 fall off is typical of p + Sj pi = 
which is a characteristic of topological defects like cosmic 
strings and textures . The inclusion of topological de- 
fects in standard cosmology has been widely considered 
(see for example 0). 

The dark radiation term in eq.JSJ now contains addi- 
tional contributions from the bulk parameters. In order 
to calculate the relative contribution of this term to the 
expansion dynamics, we write eq.© in a more convenient 
form as 



h 2 (z) = h 2 n m0 (i+z) 3 +n Aa +n d , Q (i+z) 4 -n koa0 {i+zf 

where k c g = k — a and we have defined the various en- 
ergy density components, in units of the critical density, 
in the following form. 
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The subscript denotes the value of the quantity at the 
present epoch. The BBN constrains the dark radiation 
density to be between (—1.23 and 0.11) x p r0 ||, where 
p r o is the present energy density of the photon back- 
ground with f2 r o w 2.5 x 10~ 5 . Thus the upper BBN con- 
straint imply that at the highest SN redshift z = 1.755, 
the contribution of the (positive) dark radiation term in 
eq.© is only 0.003 % of the contribution from the Sl m 
term for f2 m o as low as 0.33 (which comes from the vari- 
ous recent measurements of f2 m o as f2 m o = 0.330 ± 0.035 
at one sigma level 01 ) • Similarly for the lower constraint 
(i.e., negative dark radiation), the corresponding contri- 
bution is 0.03 %. The dark radiation term, hence, can be 
safely neglected while considering the SN la observations 
and eq.|0 reduces to 
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implying that f2 m o + Qao = 1 + ^k clt o- Note that a 
does not enter into the metric on the brane but it en- 
ters in the brane Friedmann equation eq.JSJ through the 
Israel junction conditions. This is equivalent to shifting 
the curvature index k of the standard cosmology by an 
amount —a. This term (k — ot)/S 2 can be neglected in 
the early universe. However, it dominates at the present 
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epoch and predicts significant departures from the stan- 
dard cosmology, especially the expansion dynamics of the 
universe. 

In order to test the model with the supernovae la ob- 
servations, we derive, in the following, the magnitudc- 
redshift relation. We note that the luminosity distance 
c?l of a source of redshift z, located at a radial coordinate 
distance r, is given by di, = (1 + z) So r where r can be 
calculated from the metric as 
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when k — 1 
when k = 



= sinh r when k = — 1. (12) 

The apparent magnitude to of the source is given by 

m(z) =M+5 log[P L (z)] (13) 

where M. = M — 5 logiJo+ constant, M is the abso- 
lute luminosity of the source and X>l = Hod^ is the di- 
mensionless luminosity distance. The present radius of 
the universe So, appearing in g?l, can be calculated from 
cq.jrUJ) as 
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We can now calculate, by using eas. (|llll4|l . the pre- 
dicted value of the apparent magnitude at a given redshift 
if we know the values of the parameters Q m o, S^ao, ot and 
M.. In order to test the model, we consider the data on 
the redshift and magnitude of a sample of 54 supernovae 
of type la considered by Perlmutter ct al (excluding 6 
outliers from the full sample of 60 supernovae) [ljj, to- 
gether with the recently observed supernova 1997ff at 
z = 1.755 ± 0.05 (with m co " = 25.68 ± 0.34) the highest 
redshift supernova observed so far [2(|. The best-fitting 
parameters are obtained by minimizing % 2 , which is de- 
fined by 
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refers to the effective magnitude of the «th 
supernovae which has been corrected by the lightcurve 
width-luminosity correction, galactic extinction and the 
K-correction from the differences of the R- and B-band 
filters and the dispersion <5m™ rr is the uncertainty in 
m c °". 

The model fits the data for a wide range of parameters. 
In Fig. 1, the conical volume shows the 68 % confidence 
region which gets contributions only from the closed and 
flat models. The open model contributes only at higher 
confidence levels, more than 85 %. In the case of flat 
(k = 0) model (where r2fe cff o is not necessarily zero), to, 
and hence x 2 , are independent of a which is constrained 
only through ea. (|14fl . The best-fitting parameters are 
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FIG. 1: A parameter space of interest at 68 percent confidence 
level is shown by the conical volume for a > — 1. 



obtained as f2 m o = 0.35, f^Ao = 1-12, a = 0.69 and 
M. = 23.91 with x 2 = 56.85 at 51 degrees of freedom 
(dof), which represents a fit as good as for the best-fitting 
standard model (Sl m0 = 0.87, £]ao = 1.51 and M — 23.9 
with x 2 = 56.89 at 53 dof 21] ). However, these mod- 
els are not consistent with the recent anisotropy mea- 
surements of the cosmic microwave background radiation 
(CMBR), which imply that k w However, the pa- 
rameter space is large enough (as in the case of standard 
model), as is clear from Fig 1, to accommodate with the 
CMBR predictions. Additionally, we find that the model 
also fits the data very well for the case k = for a wide 
range of parameters giving f2 m0 = 1.35, f^Ao = 1-54, 
M = 23.93 with x 2 = 58.8 (at 52 dof) as the best-fitting 
solution. Though this best fitting Q m o is higher than 
the favoured value 0.33, however the parameter space is 
sufficiently large to accommodate low fi m o values as has 
been shown in Fig 2, where we have shown the 68 % and 
95 % confidence regions for this case. Note that, unlike 
the standard model, the points in the figure are not con- 
fined only on a straight line since for a given pair f2 m o 
and f^AOj model is degenerate in a and So- For the case 
k = 0, c?l (hence to and x 2 ) is independent of a. If we 
fix f2 m0 to 0.33, then the best-fitting solution (for k = 0) 
is obtained as ft A0 = 0.7, M = 23.96 with \ 2 = 62.0 
at 53 dof. The allowed ranges of Oao at 68 % and 95 % 
confidence levels are obtained, respectively, as 0.53 — 0.84 
and 0.34 — 0.97. It should be noted that for this geomet- 
rically flat (k = 0) case, the model is not dynamically 
flat (fitotai = + ^ 1) and reduces to the standard 
flat model, Sltotai = 1 only when a — 0. 

We also notice that the exclusion of SN 1997ff from the 
sample has only a marginal effect on our results. For ex- 
ample, the global best-fitting model from the Perlmutter 
et al' sample of 54 points yields il m o = 0.34, f^Ao = 1.14, 
M = 23.91, a = 0.7 with x 2 = 56.83 at 50 dof. The 
best-fitting model for the case k = 0, from this sample, 
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yields fi m0 = 1.7, ft A0 = 1.9, M = 23.91 with x 2 = 56.9 
(at 51 dof). 
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FIG. 2: The elliptic regions at 68 % confidence level (dark 
shaded) and 95 % confidence level (light shaded) are shown 
for the geometrically flat model (k = 0). 

The present value of the Hubble parameter Hq sets 
the age of the universe. A large number of independent 
methods appear to converging on a value of Hq in the 
range (60 — 80) km s _1 Mpc -1 23], which sets the age 
of the best-fitting flat standard model in the range (11.4 
— 15.2) Gyr. If Ho is towards the lower side of this 
range, there is no serious age problem with the standard 
cosmology in the view that the age of the globular clus- 
ters t GC = 12.5 ± 1.2 Gyr Hl| and the age of Milky Way 
tuw = 12.5 ± 3 Gyr j25j. However, if Ho shifts on the 
higher side (as is claimed by the recent measurements of 
H = 72 ±7 km s^ 1 Mpc -1 by using Hubble Space Tele- 
scope 26) ), the standard model might get uncomfortable 
with its age. However, one can obtain larger age to if 
larger values of f^Ao are allowed, as is clear from Fig 3, 
where we have plotted (to) as a function of f^Ao freezing 
f2 TO o at its observed value 0.33. The expression for the 
age of the universe is the same in both models and is 
given by 
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In the favoured standard model (Qtotal — 1), can be 
increased only by increasing fl\ (i.e. by reducing fi m o), 
which does not seem likely as the recent measurements 
give very narrow range of £l m o, as mentioned earlier. 
However, there is no such constraint on the present model 
and fairly large values of Qao are allowed by the data: 
f^AO can rise as high as 0.97 at 2 sigma level (by fixing 
Jl m o = 0.33), which can increase the age as large as 14 
Gyr even for Ho as high as 72 km s _1 Mpc -1 . 

In Fig 4, we have shown the actual fitting of the best- 
fitting geometrically flat model (fixing f2 m o = 0.33) with 
the data and compared it with the favoured standard 
model. 
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FIG. 3: Age of the universe is plotted as a function of Qao for 
the cases Ho = 72 km s _1 Mpc -1 (solid curve) and Ho = 65 
km s _1 Mpc -1 (dashed curve). The value of f2 m o has been 
fixed at its observed value 0.33. The horizontal dotted line 
represents the age of the globular clusters tGC = 12.5 Gyr. 
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FIG. 4: The Hubble diagram for 55 supernovae: the theoret- 
ical curves represent the best fitting model for the case k — 
obtained by fixing f2 m o to 0.33 (solid curve) and the favoured 
standard model f2 m o = 1 — f^Ao = 0.34 (dashed curve). The 
empty circles represent the six excluded supernovae by Perl- 
mutter et al. in their primary fit C [l9|. 



We further note that the estimate of the parameter M. , 
which is usually referred to as the Hubble constant-free- 
absolute magnitude, seems to be model independent, un- 
like the other parameters (17, ) which very much depend 
on the models considered. For instance, the SN dataset, 
which we are considering here, has been fitted to a variety 
of models — the standard FRW mode ls [27| . kinemati- 
cal A- models J23 , quintessence models |2lj. auasi-steadv 
state model |2J|, etc. All give the same M. « 24, for low 
as well as high redshift supernovae j2^. This is remark- 
able and consistent with the primary assumption that 
the supernovae la are standard candles. To fix ideas, let 
us note that M. = 24 gives the corresponding absolute 
magnitude M at peak luminosity of the supernovae as 
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M = -19.3 for H = 65 km s" 1 Mpc" 1 and M = -19.1 
for H = 72 km s _1 Mpc -1 , which are in the right region 
for SN peak magnitudes ■ 

Summarizing, the modification caused by the inclusion 
of the brane scalar curvature in the action amounts to 
defining an effective dynamical curvature which is differ- 
ent from the geometric curvature. This essentially shifts 
the curvature parameter k to the effective curvature pa- 
rameter (k — a) in the Friedmann equation, hence mak- 
ing the difference even at the present epoch, which can 
be tested against the present observations. It seems to 
come out very well. This provides a welcome leverage 
which could be maneuvered for having a comfortable age 
for the universe. Our ongoing investigations indicate that 



the model seems to explain the CMBR anisot ropy mea- 
surements and radio sources data successfully |3l| . 

Although the degeneracy in the parameter space is 
large, which is due to inadequacy of the present data, the 
situation will improve when more precise data with more 
points at z > 1 become available. It would then nar- 
row down the degeneracy pinning the parameters more 
accurately. This might be accomplished by the proposed 
SNAP (SuperNova Acceleration Probe) project. 
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